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Abstract 

We derive some exact results concerning the response of strongly-interacting matter to external 
magnetic fields. Our results come from consideration of triangle anomalies in medium. First, 
we define an "axial magnetic susceptibility," then we examine its behavior in two flavor QCD 
via response theory. In the chirally restored phase, this quantity is proportional to the fermion 
chemical potential, while in the phase of broken chiral symmetry it can be related, through triangle 
anomalies, to an in- medium amplitude for vr*' ^ 27. We confirm the latter result by calculation in 
a linear sigma model, where this amplitude is already known in the literature. 
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I. INTRODUCTION 



Much attention has been focused recently on properties of matter at very high tempera- 
tures or baryon density |l| . The interest is driven by the physics of heavy-ion colhsions and 
of the core of neutron stars. Most of the discussion is focused on properties of the system 
at finite temperature T and chemical potential /i. In this paper, we are interested in the 
properties of of hot and dense matter under external magnetic field B. This question is of 
potential interest for the physics of compact objects. 

In contrast to the temperature and chemical potential, the magnetic field that can be 
achieved in nature seems to be always small compared to the strong scale (perhaps as large 
as ~ 10^^ G in magnetars This means that the effect of magnetic field on the 

medium, in most cases, can be treated as a small perturbation, and the linear response 
theory is appropriate. In this paper we derive two results related to the response of a 
strongly-interacting medium to the external magnetic field. We are interested in the axial 
current created by a uniform magnetic field B. For small magnetic fields the axial current is 
linear in B, with a proportionality coefficient which we call the axial magnetic susceptibility 
X- We show that if chiral symmetry is unbroken, then the value of x is equal to the baryon 
chemical potential, with a known numerical coefficient. This result is universal and receives 
no correction due to strong coupling. We also find that when chiral symmetry is unbroken, 
this universal result no longer holds. However in this case, we derive a relation between the 
axial magnetic susceptibility and the in-medium 7r° — > 27 amplitude in spacelike domain (a 
more precise definition is given below). In both cases the results come from the consideration 
of triangle anomalies, and hence are reminiscent to those of Refs. . The difference 



between this work and Refs. 
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13| is that here we are interested in the properties of the 
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ambient matter, not of topological defects as in Ref. 

The paper is organized as followed. Section |n] we define the quantities of interest in a 
simplified version of the real world. In Sec. IIIII we derive the exact relations. We explicitly 
check this relation in Sec. IIVI in the linear sigma model. We extend the results to the real 
world in Sec. El and conclude with Sec. IVII 
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II. BASIC DEFINITIONS 



For simplicity, we consider QCD with massless u and d quarks and neglect other quarks. 
We shall first assume nonzero baryon chemical potential, but zero isospin chemical poten- 
tial. Moreover, let us assume that electromagnetism couples to the third component of the 
isospin current ^(j'y^T^q, but not a linear combination of isospin and baryon current as in 
the real world. We shall modify our result to the real world later. We will assume no su- 
perconductivity, so the magnetic field can penetrate the matter without being confined in 
magnetic flux tubes (this may require that we work at sufficiently high temperatures). 

Let us first define the axial magnetic susceptibility. Assume we have a medium at tem- 
perature T and baryon chemical potential /i in an uniform magnetic field B. The magnetic 
field induces an axial current in the medium, 

(f) = x(T,/.)B (1) 

where, in terms of quark fields. 

That X is nonzero is permitted by symmetries. With respect to parity, both and B 
are axial vectors. Note that B and Ai has different C parity, so x must be an odd function 
of the baryon chemical potential fi. In particular, x = at zero chemical potential. In the 
model described in Sec. IIVI x is proportional to /i at small /x. 

If the medium consists of non-relativistic nucleons, then the axial current is proportional 
to the nucleon spin. The coefficient x therefore has a simple physical interpretation as 
the spin polarizability of the medium (more precisely, the the difference between proton 
and neutron spin polarizability). When nucleons are relativistic, it becomes impossible to 
separate nucleon spin from nucleon total angular momentum. However even in this case, 
the axial magnetic susceptiblity is still a well-defined concept. In fact, this axial current 
interacts, through Z° exchange, with neutrinos and modify their dispersion relation. 

We also define the in-medium coupling of the neutral pion to two photons, QnO-yy, in the 
chirally broken phase. In vacuum, the anomalous coupling of a pion to external gauge fields 
is given by a term in the chiral Lagrangian, 

_J_ .ua^Q^A^F^^^ 0=^, (3) 
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where is the gauge potential coupled to the baryon current (note that couples to the 
isospin current). The dimensionless field is normalized to have periodicity 27r. At finite 
temperature, the coupling is more subtle. As noted in Ref. 0|, there is an ambiguity with 
the zero momentum limit. We choose the following definition. Let us look at the free energy 
of a static field configuration where the 7r° field changes slowly in space, in the presence of a 
background static magnetic field B and static baryon scalar potential Aq. The free energy 
has the form 

F = S.^d,<Pf _ g^,^^_l_,vk9^^A^p,k. (4) 

Here fs is the spatial pion decay constant, and g-rrO-yy will be called the tt^ —>■ 27 amplitude. 
The free energy (0} can be thought of as arising from integrating out all degrees of freedom of 
QCD except the Goldstone boson, and restricting to the lowest Matsubara frequency u = 0. 
At zero temperture fs = fn and g-n-o^^y = 1, but in general both and S'ttO^^ are functios of 
temperature and baryon chemical potential. We also know that ^ at the second-order 



chiral phase transition, with the critical exponent of Josephson's scaling 



Q 



. The way we 



define 5'7r077 corresponds to the — > 27 amplitude in the spacelike region of Ref. 
To summarize, our results are 
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• When the point (T, /i) lies in the chirally restored phase, x is directly proportional to 
the chemical potential. 

The numerical coefficient l/(47r^) is exact and is related to triangle anomaly. 

• When chiral symmetry is sponteneous broken, the relation between x and the anomaly 
is lost. However, there is an exact equation relating the susceptibility and the in- 
medium amplitude of 7r° 27. Namely, 

4vr^^ + g^rO-y^iT, /i) = 1. (6) 

Here g^^o^^iT) is the vr*^ 27 amplitude defined above. 

We note that the first part of our results, which concerns the chirally restored phase, has 
been observed in Ref. 0]. In addition, it is only a slight variation on the result that is found 
in Ref. 7[ , where it was determined that the magnetic susceptibility of the electric current is 
proportional to a chemical potential for fermion chirality. This relation can also be checked 



explicitly in models of strongly-interacting field theory with gravity dual description j8|. 
However, to our knowledge, the case with spontaneous breaking of chiral symmetry has 
never been considered before; hence the second part of our results is new. We first show 
the validity of these relations in a general setting. Then we shall verify them explicitly in a 
model with anomaly, namely the linear sigma model. 



III. EXACT RELATIONS 

To derive the exact relations, we consider a three-point correlation function of the axial 
current j^, the isospin current j^, and the baryon current B^, 



where j^^ is defined in Eq. (j2I) and other currents are defined as follows 

f = 1 (uji^u - d-f^d) , 3" = ^ (m7'^m + dj^d) . (8) 

We can always define the correlator so that the triangle anomaly resides entirely in the 
derivative of the axial current: 

p^T^-\p,q) = -^e^'^^p.qp, (9) 
q^V^^'ip, q) = (px + qx)r''''\p, q) = 0. (10) 

In this paper we will be interested only in static (time-independent) problems, therefore 
we shall set Pa = qo = 0. Moreover, the baryon chemical potential couples to the zeroth 
component of B\. Thus the quantity of interest for us will be 

zr^°(p,q) = y"rfS Ae-*P-™y(A^(a;)l^^(?/)5°(0)). (11) 

On the other hand, the axial magnetic sucsceptibility x is related to the low-momentum 
behavior of a two-point function. Indeed, the axial current created by a background elec- 
tromagnetic field is 

{3''{x)) = -ij d'^yG^.n^ - y)My), (12) 

where 

G^';{^-y) = {jn^)fiy))- (i3) 
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In order to reproduce Eq. (^1, the infrared behavior of the G51 correlator must be as follows: 



d'xe-^^-'^ifWiO)), = xe'^'Pj + 0{p'). (14) 

In Eq. (fT^ we emphasize that the average is taken at nonzero chemical potential /i. 
Differentiating Eq. (fT^ with respect to /i we get 

r'=0(p, _p) = _M^e*jy + o(/). (15) 

Let us now look at the structure of the correlator rjjo(p, q) in the regime of small p and 
q. It changes sign under parity. This static correlator is not singular in the chirally restored 
phase, and may have a pion pole in the chirally broken phase. In light of this, the general 
form of the three point function is, 

r^^°(p, q) = Cie'^'^q" + Cs^e'^'^V- (16) 

Since F is even under C, the dimensionless constants Ci and C2 must both be even under 
C parity, with C2 vanishing in the chirally restored phase. 
From the relation of triangle anomaly we find 

C, + C. = ^. (17) 

Consider first the case when chiral symmetry is restored. Then as the singular term in 
Eq. is absent and Ci = l/(47r^). But by comparing Eq. with Eq. (jl5j) . we find 

Requiring x to be an odd function of fi, we determine 

xif^^T) = £-^, (19) 

which is the first part of our result. 

Now we turn to the chirally broken phase. The singular term in Eq. (fT^ comes from 
the Feynman diagram with an intermediate pion line, which can be computed using the free 
energy as an effective Lagrangian. This leads us to 

C2 = (20) 



Equation ()17|) then implies 
which is the second part of our result 



4vr'^ + ^^o^^ = 1, (21) 



IV. EXAMPLE: LINEAR SIGMA MODEL 



Since QCD is strongly coupled for temperatures below the chiral phase transition, we 
cannot directly check the formula (jH} there (although it should be possible to verify it 
in the high-density phase of three-flavor QCD, where chiral symmetry is broken at weak 
coupling) . We shall instead verify this formula in a weakly coupled field theory with anomaly, 
namely the linear sigma model. This model was employed before to understand the effects 
of temperature on anomaly j9j. Furthermore, the finite-temperature vr'' 27 amplitude 
has been computed in this model in various kinetimatic limits, including the limit where 
the outgoing photons are at zero frequency . Thus, we can confirm our result (jHI) by 

calculating the axial magnetic susceptibility in this model. 

The model is given by the Lagrangian 

The couplings g and A are small. The expectation value of a is f, which is temperature 
dependent (being equal to a/ /i^/ A in vacuum.) The covariant derivative, = 9^ + iqA^j,, 
is that of f/(l)EM- In the phase with chiral symmetry breaking, fermions ("constituent 
quarks") have mass m = gv. We shall calculate the axial current induced by a baryon 
chemical potential /i^ on the background of a constant magnetic field, using the single- 
particle Hamiltonian in the regime where T ^ /x, m, VeB. 

Placing a static, homogeneous magnetic field pointing in the z direction can be accom- 
plished by means of the vector potential. A'^ = (0, 0, Bx, 0). The fermion spectrum can be 
found by solving the Dirac equation ■ DQ = EQ with D'^ = — iqA'^, where, q is the 
electric charge of the quark Q = {u, d), = —qd = |- The axial- vector current can then be 
directly calculated as a thermal expectation value of (j21). 

Within this setup it is convenient to parametrize the quark wave functions as, 

f M^) \ 

(j)l{x) 

\ M^) J 

Because of our choice of coordiante system, "p^" will always appear in the same component 
of the Dirac equation as "ei?a;" in the linear combination py + eBx. Thus, by making a 



Q{x,y,z) 



^i{pyy+p.z)^ (23) 



convenient coordinate shift, x — > x — || we can eliminate the exphcit appearance of "p^" in 
the Dirac equation. We have four first order coupled DEs, from which we can obtain two 
second order equations for ■02 and 02 separately: 



[V' - {qBxf + -pl-m^ + qB]^^) = 0, 



'.{x) = 0. 



(24) 
(25) 



The other components can be expressed via ■02 and 02, as 

{ E - p^)^2 - ^02 

E"^ — pI — IV? ' 
{E + p^)(f)2 - rmp2 



^ i\S + qBx]- 
01 = -i[V + qBx 



(26) 
(27) 



E'^ — pI — vn? 

The components ^02 (a:;) and 02 (x) are clearly the eigenfunctions of a harmonic oscillator, 
while the other components are obtained by acting a lowering operator on linear com- 
binations of ^02 and 02- The energy eigenvalues are thus the familiar Landau levels, 
E — ±-y/p^ + vn? + 2nqB. The constituent quark wave function, properly normalized, is 



/ ^2(E±p^nn-i{VqBx)\ 



Q{x,y,z) 



1 



aHn{\/qBx) 

^ '^'"-'>(^-p^) n,,_,i^x) 



-qBx/2 i{pyy+pzz) 



(28) 



\ PHniVqBx) 
The Landau levels, n, replace the Px quantum number. 

The H are the Hermite polynomials, with the addition that, forn = we define H-i = 
For each n > there are two eigenstates with each choice of the sign of E: 

( m — \/1nqB \ 

qB ^ 




167r 



VE{E-p, 
E-p, 



\ 



fqB_ 
V167r 



1/4 



E 



m + ^/2nqB 
V ^E{E + p,) ) 



(29) 



E ) 

Meanwhile, for n = there is only one positive energy eigenstate and one negative energy 
eigenstate. 




qBy 
An 



( lE + p, \ 
E 



E-p, 



(30) 



\y E / 
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Now we calculate, 

{AD = {Q^.^'Q)^,T - {Qi'i'Q)o,T. 
by filling up the fermion energy levels with the Fermi-Dirac distribution function, 

1 1 , 



n,A sgn{E) ' 



m' + \4>iV-\m[ 



z. 



(31) 



(32) 



Here /iq = |/i is the quark chemical potential, and the prefactor (gi?)/(27r) arises from the 
degeneracy of the energy eigenstates in py. The sum is taken over all Landau levels n, both 
sign of energy E and all polarization A. For regularization we also subtracted the value of 
the axial current at /i — 0, which is zero by C parity. For n > the sum over polarizations 
read 



-2m^}, (33) 



which is an odd function of and contributes nothing to the axial current after integration 
over dpz- It is very easy to see that the n > Landau levels do not contribute to the 
axial current in two limits, when the fcrmions are massless and when the fermions are 
nonrelativistic (very massive) , so it it not entirely surprising that they do not contribute for 
any value of m. 

The contribution from the lowest Landau level is considerably simpler, since here 



[|V'l|-|V'2|'+|0l|'- 10211 = 



= 1. 



(34) 



All that remain is to calculate the integral over the statistical factor. This can be done 
analytically at large temperatures by expanding to first order in the small quantities, /i/T 
and /T'^. For the sake of brevity we use n{^) = (e^ + 1)~^, and n'(^) = d^n{^), with 
^ = Pz/T. In particular. 



1 



d 



-n{(3E) = -I3nq 



n'(0 + ^^n"(0 



(35) 



The sum over flavors yields — = 1 in place of "g" in the prefactor qB /{27r). Also, a 
factor of A^c = 3 is gained from the sum over colors. Thus we find 



4^ 



m 



rp2 



(36) 
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Performing the integral as in Ref. we obtain 



dx 



Note that we recover Eq. (0) for m = 0. 

On the other hand, the result of Refs. 0,0, in our languague, corresponds to 

7C(3)m2 



(37) 



We see that 



in accordance with Eq. (jHl). 



77 



47r2T2 



4^ ^ + 9^° 



77 



(38) 



(39) 



V. ISOSPIN CHEMICAL POTENTIAL, REAL- WORLD EM COUPLING 

In general, the system could contain an isopsin chemical potential as well. This section 
expresses the results of the last two sections, as they partain to two flavor QCD with non- 
zero /i/ and /is, in an applied magnetic field coupling properly to f/(l)EM- The f/(l)EM 
vector current can be expressed in the flavor basis, as a linear combination of the isospin 
and baryonic generators. Specifically, Jq = j^j^ = e(|jg + jf). We now expect, 

(f) = x(T,/xi,/iB)BQ. (40) 

We will see that x{T, f^i, fis) is a linear combination of /ij and /xb- 

In the chirally symmetric phase, we can follow section ITTTl but replacing the A" coupling 
to isospin with one coupling to electric charge, Aq. Then the response of the chiral current 
can be expressed in terms of the correlator, Gsg as 



{f^{x)) = -zJ d'yG^'Q{x - y)A,Q{y). (41) 

For this system, we must have 

j d^x e'^-'^GfQix) = xe'^'^Pj + 0{p^). (42) 

Now, two different anomaly relations - one for (j^^jgE^) and one for {j^^igV^) - give 
respectively, 

^ = — ^ = ±- (43) 

dpiB 47r2' dfjii 87r2' ^ ' 
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For the phase with broken chiral symmetry, the analysis of section IIIII still holds as well, 
but with the distinction that there would be two separate diagrams with propagating pions, 
necessitating the introduction of two amplitudes, g-^o^j and g.„o^B defined analogously with 
Qt^o^^. Again, restricting the anomaly relation to different space time components results in 
two different equations. 



871 h g^o^j = 1, Att h g^o^B = 1- (44) 

C/i/ OfJ'B 

The 7r° 77 amplitude can be defined as a linear combination of gTrO^j and g-^o^B'- 



47r2 2 47r2 avr 



{gnO-yy IS normalized to 1 at zero temperature.) From Eqs. ()44|1 we find 

4^'fl^ + ^I^V^?-°77 = l- (46) 



This result could again be confirmed for the sigma model by noting that the sum over 
flavors should be performed with fi^ 7^ f^d, resulting in "g" in the prefactor being replaced 
by qufJ'u - qdfJ'd- Envoking hb = + yUd) and hj = - n^) we find 



gn/in - QdfJ-d = 2 yj'B + ^1^1 j ■ (47) 

Thus, 

„. ^ ji, . 1„) [^-m'^,0{!^).0(^).0 (^)] . ,48, 
From the results of Refs. QQ, gnO^^ = 7C{3)m'^/{A7c^T^), one can verify Eq. 



VI. CONCLUSIONS 



We have shown that there is a close connection between the response of strongly interact- 
ing matter on external magnetic field and the axial anomaly. By considering the properties 
of the three point function of isovector, axial isovector, and baryon currents in the pres- 
ence of non-vanishing baryon chemical potential and QED magnetic field we were able to 
show that the axial- magnetic susceptibility, x(/^) T) is directly proportional to the baryon 
chemical potential in the absence of Goldstone modes carrying the axial current. Alekseyev 



11 



et. al. previously found the same to be true for massless QED, using quite general meth- 
ods 0. This result follows from the fact that the anomaly coefficient, which determines 
the constants value, is not renormalized, and receives no finite-temperature contibution. In 
the presence of massless pions, this direct relation no longer holds but one still can relate 
T) to the anomaly coefficient through the amplitude for tt*^ — > 77. We confirmed the 
second relation for a particular case of weakly-coupled linear sigma model. 

These results may be applicable to the study of compact objects such as neutron stars, 
where both baryonic chemical potential, and magnetic field may be large. Specifically, the 
self energy of neutrinos is affected by interaction with the axial isovector current through 
exchange Q- Calculation of 5'^o^^(/i,T) in nuclear matter is complicated by the need 
to deal with singularities arising from particle-hole interactions, but our results could be 
easily employed to find this neutrino self energy contribution in deep cores of neutron stars, 
if chiral symmetry is restored there. In a strong magnetic field the momentum distribution 
of neutrinos is asymmetric already in equilibrium, so they will stream out in asymmetric 
fashion, giving rise to a small contribution to pulsar velocities. The presence of an axial 
current in matter will not affect oscillations between active neutrinos, but does change the 
oscillations between an active neutrino and a sterile one. 

We end the paper by noting that currently we lack an understanding of the critical 
behavior of Qj^o^^ near the second-order chiral phase transition. While it is natural that this 
coefficient goes to zero smoothly at the phase transition, the question about the critical 
exponent remains open. 

This work is supported, in part, by DOE grant DE-FG02-00ER41132. D.T.S. is sup- 
ported, in part, by the Alfred P. Sloan Foundation. We thank A. R Zhitnitsky for comments 
on an earlier draft of this manuscript. 
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